NASA 

Technical  Memorandum  105653 


AVSCOM 
Technical  Report  91-C-027 


AD-A252  380 


Optimum  Design  ot  a  uearbox 
for  Low  Vibration 


DTIC 

ELECTE 

JUL021992 

A 


Katsumi  Inoue,  Dennis  P.  Townsend, 
and  John  J.  Coy 
Lewis  Research  Center 
Cleveland,  Ohio 


I  This  document  has  been  approved 
[  lot  public  release  and  solo;  its 
I  dijljibution  is  unliinitsd. 


92-17101 

iiai  lait 


Prepared  for  the 

Sixth  International  Power  Transmission  and  Gearing  Conference 
sponsored  by  the  American  Society  of  Mechanical  Engineers 
Scottsdale,  Arizona,  September  13-16,  1992 

92  b 

f\l/NSA 


SYSTEMS  COMMAND 


E-G622 


OPTIMUM  DESIGN  OF  A  GEARBOX  FOR  LOW  VIBRATION 

Katsumi  Inoue^ 

Tohoku  University 
Sendai,  Japan 

and 


Dennis  P.  Townsend  and  John  J.  Coy 
National  Aeronautics  and  Space  Administration 
Lewis  Research  Center 
Cleveland.  Ohio  44135 


Accesion  For 

NTtS  CRA&I 
one  TAB 
U..a:M.ouncec; 
Ju;ttfication 


By . . . 

Di,  1  ibuiion/ 


Availabilll'/ 


J 


ABSTRACT 

A  computer  program  was  developed  for  designing  a  low  vibra¬ 
tion  gearbox.  The  code  is  based  on  a  finite  element  shell  analysis 
method,  a  modal  analysis  method,  and  a  structural  optimization 
method.  In  the  finite  element  analysis,  a  triangular  shell  element 
with  18  degrees-of-freedom  is  used.  In  the  optimization  method,  the 
overall  vibration  energy  of  the  gearbox  is  used  as  the  objective 
function  and  is  minimized  at  the  exciting  frequency  by  varying  the 
finite  element  thickness.  Modal  analysis  is  used  to  derive  the 
sensitivity  of  the  vibration  energy  with  respect  to  the  design 
variable.  The  sensitivity  is  representative  of  both  eigenvalues  and 
eigenvectors.  The  optimum  value  is  computed  by  the  gradient 
projection  method  and  a  unidimensional  search  procedure  under  the 
constraint  condition  of  constant  weight. 

The  computer  code  is  applied  to  a  design  problem  derived  from 
an  experimental  gearbox  in  use  at  the  NASA  Lewis  Research 
Center.  The  top  plate  and  two  side  plates  of  the  gearbox  are  rede¬ 
signed  and  the  contribution  of  each  surface  to  the  total  vibration  is 
determined.  Results  show  that  optimization  of  the  top  plate  alone  is 
effective  in  reducing  total  gearbox  vibration. 

INTRODUCTION 

The  vibration  of  a  power  transmission  gear  system  is  initiated 
by  the  gear  mesh.  The  vibration  energy  is  transmitted  to  the  gear¬ 
box  housing  through  the  shafts  and  bearings,  and  the  gearbox 
housing  radiates  the  structure-borne  noise.  The  vibration  and  gen¬ 
erated  noise  are  influenced  by  manufacturing  precision,  alignment 
errors  of  the  gears,  and  the  design  of  the  gearbox.  These  factors 
affect  the  total  quality  of  the  system  and  are  becoming  more 
important. 

In  order  to  estimate  the  vibration  of  a  gear  system,  it  is  neces¬ 
sary  to  understand  the  excitation  characteristics  and  to  evaluate  the 
dynamic  behavior  of  the  system.  With  respect  to  excitation  char¬ 
acteristics,  Kubo,  Kiyono,  and  Fujino  (1986)  proposed  the  concept  of 
the  total  vibration  exciting  force  that  integrates  the  effects  of 


manufacturing  and  alignment  errors  and  periodical  change  of  tooth 
mesh  stiffness.  However,  the  evaluation  of  the  system’s  dynamic 
behavior  is  related  to  the  modeling  of  the  system.  Some  approaches 
presented  recently  include  the  analysis  of  multistage  gear  systems 
coupled  with  gearbox  vibrations  (Choy  et  al.,  1991)  and  the  analysis 
of  vibration  transmission  in  a  single-stage  gearbox  by  the  building 
block  approach  (Takatsu  et  al.,  1991).  These  are  not  simple  simu¬ 
lations.  However,  once  the  vibration  of  a  gearbox  is  analyzed,  the 
sound  radiation  from  the  gearbox  can  be  estimated  by  using  the 
results  as  input  data  to  BEMAP  (boundary  element  method  for 
acoustic  prediction)  (Seybert  and  Wu,  1989)  or  a  similar  boundary- 
element-based,  acoustic-intensity  analysis.  Consequently,  the 
vibration  and  noise  of  a  gearbox  can  be  estimated  to  some  extent. 

Generally,  the  vibration  of  the  gear  system  is  sensitive  to  the 
dynamic  characteristics  of  the  gearbox.  The  reduction  of  noise  levels 
is  linked  to  reduction  of  the  gearbox  vibration  levels.  In  contrast  to 
progress  made  in  the  dynamic  analysis  of  gear  systems  mentioned 
previously,  the  systematic  procedure  for  reducing  gearbox  vibration 
levels  is  not  adequate.  The  technique  to  avoid  a  resonance  by  shift¬ 
ing  the  natural  frequency  is  often  used  to  reduce  vibration.  A  mass 
or  stiffener  is  added  at  an  efficient  position  on  the  gearbox.  Al¬ 
though  this  can  be  an  effective  method,  it  is  not  easy  to  determine 
the  best  position  on  the  gearbox  or  to  estimate  a  suitable  natural 
frequency  shift.  From  this  viewpoint,  the  authors  have  proposed  a 
method  to  design  plates  (Inoue,  Kato,  and  Ohnuku,  1990)  and  thin- 
plate  structures  (Inoue,  Townsend,  and  Coy,  1992).  The  overall 
vibration  energy  is  adopted  as  the  parameter  for  estimating  the 
quality  of  a  design  and  is  minimized  at  a  given  frequency  by  varying 
the  finite  element  thickness.  The  optimization  problem  is  solved 
under  the  constraint  of  constant  weight.  In  this  paper  the  method  is 
applied  to  the  design  of  a  gearbox  made  from  steel  plates,  and  the 
vibration  reduction  is  demonstrated  analytically.  The  exciting  forces 
applied  at  the  bearing  locations  are  assumed  to  be  known,  although 
influenced  by  the  gearbox’s  dynamic  characteristics  determined  in 
the  design  process.  The  purpose  of  this  research  was  to  design  a 
low-vibration  gearbox. 
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OPTIMUM  DESIGN  PROCEDURE 

Minimization  of  Overall  Vibration  Energy 

The  theoretical  basis  of  the  proposed  method  is  given  by  Inone, 
Townsend,  and  Coy  (1992)  and  is  summarized  here.  The  equation 
of  motion  of  a  structure  is  given  in  the  matrix  form  as  follows; 


IMl  {u(t)}  +  !C;  {u(t)}  +  IK)  {u(t)}  =  {f(t)}  (1) 

where  (M)  is  the  mass  matrix,  (C)  is  the  damping  matrix,  [K!  is  the 
stiffness  matrix,  {u(t)}  is  the  time-varying  displacement  vector,  and 
{f(t)}  is  the  time-varying  exciting  force  vector.  Proportional  damp¬ 
ing  is  assumed,  and  the  damping  matrix  is  obtained  by  the  linear 
combination  of  the  mass  matrix  and  the  stiffness  matrix: 

[C]  =  ajM]  djK]  (2) 

where  o  and  ^  are  the  coefficients  to  define  the  proportional 
damping  matrix.  The  modal  matrix  [$  j  is  defined  by  the  assembly 
of  normalized  eigenvectors 
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where  R  and  1  are  the  real  and  imaginary  parts  and  {...}^  de¬ 
notes  transpose  matrix.  The  vibration  energy  T  of  the  structure  is 
given  by 


T  1  {u  ’FiM]  {u} 

2 

=  -  {'i  :M1  ;$1  {«} 

2 

=  Mrr{^} 
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[$!  =  {d,}  ,  .  .  .  ,  {<*4. 


where  {U  *}  and  {{  ’}  show  the  conjugate  complex  of  {U}  and  {f}. 
The  energy  is,  therefore,  represented  by  the  following  expression: 


where  n  indicates  the  number  of  eigenvectors  considered.  When 
the  structure  is  excited  by  a  harmonic  force  of  angular  frequency  w, 
the  force  vector  and  vibration  displacement  vector  are  given  by 


(10) 


{f(t)}  =  {F}ei“‘ 
{u(t)}  =  (U}ei“‘  = 


(4) 


where  {F}  is  amplitude  of  exciting  force,  (U)  is  amplitude  of  dis¬ 
placement,  (f)  is  the  vector  defined  as  {U}  =  and  j  is  the 

imaginary  unit.  Substituting  Eq.  (4)  in  Elq.  (1),  and  premultiplying 
it  by  [$|^,  the  n  independent  algebraic  equations  are  obtained  to 
evaluate  the  vector  {f}.  The  rth  component  is  given  in  the  fol¬ 
lowing  form  by  solving  the  equation 


^r 


/  ((■'r  -  +  j“(“  t-  ^-'r))  =  1.  2, 


(5) 


In  the  expression,  A^.  denotes  the  rth  eigenvalue.  Therefore,  the 
amplitude  {U}  for  the  vibration  displacement  is  obtained  by  sub¬ 
stituting  ({}  into  Eq.  (4). 

The  vibration  velocity  {u(t)}  is  given  by  the  similar  form  as 
shown  in  Eq.  (4). 


The  sensitivity  of  the  vibration  energy  with  respect  to  the  design 
variable  Xj  (i  =  1,  .  .  .  ,  m)  is  given  by 


iT  =  r 


(H) 


Since  a,  f,  and  (F)  in  Eq.  (8)  are  independent  of  the  design 
variable,  the  sensitivity  of  the  vibration  energy  is  represented  by  a 
function  of  the  sensitivities  of  the  eigenvalue  and  the  eigenvector. 
The  evaluation  of  the  sensitivities  of  both  eigenvalue  and  eigenvector 
is  based  on  the  method  proposed  by  Fox  and  Kapoor  (1968). 

The  optimum  solution  is  obtained  by  minimizing  the  vibration 
energy  under  the  constraint  of  constant  weight.  The  element  thick¬ 
ness  is  adopted  as  the  design  variable  and  limited  by  specified  mini¬ 
mum  and  maximum  values.  The  gradient  projection  method 
(Rosen,  1961)  weis  used  for  the  optimization  procedure  in  this 
research.  The  value  of  the  energy  at  every  iteration  is  saved  and 
used  to  determine  convergence.  The  convergence  criteria  is  1/1000 
relative  variation. 


(u(t)}  =  {u}el“‘ 


(6)  Optimum  Design  Program 


The  velocity  can  al.so  be  derived  directly  by  differentiating  the 
displacement  in  Eq.  (4);  therefore,  the  following  equation  is  ob¬ 
tained  by  equating  the  differentiation  to  Eq.  (6): 

If}  =  jw{f} 


A  computer  program  was  developed  based  on  the  procedures 
previously  described.  The  program  consists  of  two  main  pans.  The 
first  part  is  the  analyzer,  which  is  composed  of  the  finite  element 
analysis,  the  eigenvalue  analysis,  and  the  evaluation  of  the  vibration 
energy  using  modal  analysis.  The  second  part  is  the  optimizer, 
which  includes  the  sensitivity  analysis  and  search  procedure  for  the 
optimum  value. 


Therefore,  the  real  and  imaginary  parts  of  the  rth  component 
are  obtained  as  follows: 


In  the  finite  element  analysis,  the  program  is  coded  to  make  the 
design  program  as  compact  and  flexible  as  possible.  If  the  total 
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degrees-of-freedom  is  constant,  using  a  small  number  of  high- 
precision  elements  generally  leads  to  good  results.  However,  the 
acceptable  degree-of-freedom  of  the  element  should  be  within  limits 
in  this  research  because  several  elements  are  necessary  to  model  the 
three-dimensional  complex  surface.  Therefore,  a  triangular  shell 
element  with  18  degrees-of-freedom  is  used.  The  triangular  shell 
element  is  formed  by  combining  the  elements  for  plate  bending  and 
plane  problem  analyses  given  by  Zienkiewicz  and  Cheung  (1967). 

The  consistent  mass  matrix  used  for  this  analysis  is  evaluated  using 
the  numerical  integration  method.  The  accuracy  of  the  finite 
element  analysis  is  briefly  discussed  in  the  Appendix. 

The  program  is  configured  on  a  Digital  Equipment  Corporation 
VAX  cluster  at  NASA  Lewis  Research  Center  and  executed  inter¬ 
actively  for  small  models.  For  large  models  the  program  is  run  on 
the  Cray  XMP/28  supercomputer.  The  program  is  written  in 
FORTRAN  77  with  no  special  software  or  hardware  requirements. 

MODELING  AND  EIGENVALUE  ANALYSIS  OF  GEARBOX 

Figure  1  shows  the  gearbox  used  as  a  design  model.  This 
gearbox  is  used  in  the  NASA  Lewis  gear-noise  test  facility  (Oswald 
et  al.,  1991)  and  is  made  of  steel  plates.  The  origin  of  a  Cartesian 
coordinate  frame  is  placed  at  the  center  of  the  bottom  plate,  and  the 
positive  Z  axis  is  defined  upward  and  normal  to  the  bottom  plate. 
The  Y  axis  is  parallel  to  the  direction  of  the  shafts.  The  overall 
dimensions  of  the  gearbox  are  approximately  324  by  248  by  270  mm 
(dimensions  in  X,  Y,  and  Z  directions).  The  four  side  plates  are 
welded  together  and  to  the  bottom  plate.  The  Lop  plate  is  bolted  to 
the  side  plates.  To  simplify  the  formulation,  the  gearbox  is  modeled 
as  a  simple  box  by  using  256  triangular  shell  elements  with  all 
surfaces  continuously  connected.  The  stiffeners  and  the  four  circular 
holes  for  the  bearings  are  neglected  in  the  model  as  is  the  supportive 
effect  of  the  shafts.  These  simplifications  may  be  reexamined  and 
more  elements  used  for  a  precise  analysis,  while  still  conserving  the 
basic  concept  of  the  proposed  method.  The  number  of  elements  is 
fixed  in  this  research  only  to  limit  computer  time.  The  four  corners 
of  the  bottom  are  fixed.  Figure  2  shows  the  mesh  pattern. 


Figure  2. — Mesh  pattern  used  for  the  optimum  design  of  NASA 
gearbox. 


The  six  surfaces  are  identified  by  the  numbers  I  to  6  as  follows: 
surface  1,  the  top;  surfaces  2  and  4,  the  plates  at  Y  "  124  and 
-124;  surfaces  3  and  5,  the  plates  at  X  =  162  and  -162;  and  sur¬ 
face  6,  the  bottom.  The  thickness  of  the  portions  for  the  bearing 
housings  at  -121  ^  X  <  82  and  66.5  <  Z  <  193.5,  is  22.4  mm;  the 
other  elements  arc  6.3  mm  thick.  The  constants  used  in  this  re¬ 
search  are  modulus  of  elasticity  -  206  GPa,  Poisson’s  ratio  =  0.3, 
mass  density  -  8000  kg/m^,  coefficients  for  composing  proportional 
damping  matrix  o'  =  1.0  sec  ',  and  ?  —  5.0  <  10'^  sec. 

The  eigenvalue  problem  is  solved  first  to  obtain  the  dynamic 
characteristic  of  the  gearbox.  The  lower  10  natural  frequencies  and 
mode  shapes  are  summarized  in  Table  1.  The  mode  shape  of  each 
surface  is  represented  as,  for  example,  (m„,  m^),  which  indicates  that 
m,i-land  nu  -1  nodal  lines  arc  observed  in  X  and  Y  directions. 
The  calculated  fundamental  natural  frequency  of  475  Hz  is  approxi¬ 
mately  4.4  percent  less  than  the  result  of  measurement  in  the  report 
by  Lim,  Singh,  and  Zakrajsek  (1989).  Some  other  measurement 
data  are  shown  in  this  report,  but  there  are  no  modes  that  coincide 
with  the  calculated  results.  This  may  be  caused  by  neglecting  the 
effects  of  stiffeners  and  shafts. 


VIBRATION  OF  GEARBOX  RESULTING  FROM  UNIT 
HARMONIC  FORCES 

Unit  harmonic  forces  IN  are  applied  to  the  four  positions  of  the 
bearings  in  the  direction  of  the  line  of  action,  and  the  total  vibration 
energy  is  calculated  for  the  exciting  frequency,  which  is  given  at 
every  10  Hz  between  500  to  1500  Hz.  The  energy  is  evaluated  on 
the  basis  of  the  lower  30  modes.  Since  the  natural  frequency  of  the 
30th  mode  is  2448  Hz,  which  is  higher  than  the  considered  frequency 
range,  30  modes  is  enough  for  this  calculation.  The  vibration  energy 
is  shown  by  the  broken  line  in  Fig.  5.  Ten  peaks  are  observed  in  the 
range  al  540,  580,  610,  770,  940,  1240,  1280,  1320,  1400,  and 
1460  Hz,  respectively.  Since  the  energy  is  calculated  at  every  10  Hz, 
the  peaks  do  not  coincide  exactly  with  the  natural  frequencies,  some 
of  which  are  shown  in  Table  1.  However,  they  indicate  the 
resonance  approximately. 
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TABLE  l.-LQWER  to  natural  FREQUENCIES  AND  MODE  SHAPES  OF  GEARBOX 


Nttural 

frequency 

duotber. 

Fre- 

queacy 

lU. 

lU 

Surface  number,  mode  ithape 

i(“'..'»,i 

2(m^,m,) 

3(my.mJ 

4(m„,m,} 

5(«>,.“.) 

1 

473 

‘(1.1) 

‘(M) 

•(1.1) 

‘d.n 

‘(1,1) 

•(1.1) 

2 

540 

(1.2) 

‘(1.1) 

(2.1) 

•d.i) 

(2.1) 

(1.2) 

3 

543 

Ni.i) 

(1.2) 

(1.2) 

(1.2) 

(1.2) 

•(1.1) 

4 

582 

‘(1.1) 

•(1.1) 

‘(1.1) 

■(1.1) 

‘(1.1) 

•(1.1) 

5 

607 

(2.1) 

(2.1) 

•(l.l) 

(2,1) 

‘(1,1) 

(2.1) 

6 

768 

•(1.1) 

•(1.1) 

•d.i) 

•d.i) 

•(1.1) 

*(1.1) 

7 

942 

(2.1) 

‘(1.2) 

(1.2) 

‘(1.2) 

(1,2) 

(2.1) 

8 

994 

(2.1) 

(2.1) 

•(1.1) 

(2,1) 

‘(1.1) 

(1.1) 

9 

1041 

(1,2) 

‘(1.2) 

(2.1) 

‘(1,2) 

(2,1) 

(1.2) 

10 

1101 

(l.'O 

‘(1.2) 

(2,1) 

■(1.2) 

(2.1) 

(1.2) 

*C«>nve*  oul««rd. 
^CoBcavc  outwtrd. 
'Oblique  nod^l  line 


Generally,  the  product  of  an  exciting  force  and  the  correspond¬ 
ing  displacement  gives  the  work  done  by  the  force.  Since  the  work 
done  represents  the  energy  input  into  the  vibrating  structure,  it 
must  have  a  strong  correlation  with  the  vibration  energy  defined  in 
this  research.  This  correlation  is  demonstrated  approximately  for 
the  aforementioned  frequencies  of  the  gearbox  in  the  following 
manner. 


The  mode  shape  obtained  by  the  eigenvalue  analysis  is  not  the 
actual  displacement,  but  the  relative  displacement  at  every  node. 
However,  all  the  eigenvectors  calculated  in  this  research  are  normal¬ 
ized  with  respect  to  the  mass  matrix.  Therefore,  if  every  eigenvector 
component  is  listed  at  a  certain  node,  it  might  represent  the  relative 
linear  and  rotational  components  of  the  displacement  at  the  node. 
The  components  at  the  position  of  the  bearings  are  selected  from  the 
results  of  the  calculation  and  normalized  with  respect  to  the  maxi¬ 
mum  value  to  obtain  the  displacements.  Figure  3  illustrates  the  dis¬ 
placements  at  bearing  1.  The  abscissa  indicates  the  mode  order 
from  1st  to  30th.  For  example  (Figs.  3(a)  and  (c)),  the  maximum 
displacements  in  the  X  and  Z  directions  originate  in  the  27th  and 
28th  mode  shape,  respectively. 

The  product  of  the  X  and  Z  components  of  the  applied  unit 
load  and  the  corresponding  displacements  at  the  bearing  positions  is 
defined  here  as  the  work  done  by.  the  load, 

Li  =  iu,  -  Uji  sin  +  |w,  -  Wjl  cos  (12) 


where  u  and  w  arc  the  displacements  in  X  and  Z  directions; 

is  the  operating  pressure  angle;  subscripts  1  and  2  indicate  the 
position  of  bearings  1  and  2,  respectively.  The  displacements  Uj 
and  Wj  are  taken  from  the  results  at  the  aforementioned  10  fre¬ 
quencies:  the  2nd,  4th,  5th,  6th,  7th,  11th,  12th,  13th,  15th,  and 
17th  modes  in  Fig.  3.  A  similar  method  is  used  to  obtain  Uj  and 
Wj,  and  the  work  done  is  calculated  from  Eq.  (12).  The  work  done 
is  not  the  actual,  but  the  relative  work.  Figure  4  shows  the  com¬ 
parison  of  the  work  and  the  calculated  vibration  energy  of  the  gear¬ 
box.  The  energy  is  approximately  proportional  to  the  work.  This  is 
as  expected  and  suggests  that  the  relative  vibration  magnitude  can 
be  estimated  approximately  from  the  work. 


OPTIMUM  DESIGN  OF  GEARBOX 
Design  of  the  Top  Plate 
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The  vibration  energy  of  the  gearbox  caused  by  unit  exciting 
force  IN  with  initial  thickness  is  shown  by  the  broken  line  in  Fig.  5. 
The  abscissa  indicates  the  exciting  frequency,  and  the  peaks  repre¬ 
sent  the  resonance.  In  the  first  example,  the  top  plate  of  the  gear¬ 
box  is  designed  through  the  energy  minimization  at  1260  Hz.  This 
frequency  coincides  with  the  gear  mesh  frequency  for  a  28-tooth  gear 


(c)  Displacement  along  Z  axis  and 
rotation  about  X  axis. 

Figure  3. — Relative  displacements  at  the 
positions  of  bearing  1  (normalixed  by  the 
maximum  displacement). 
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Work  done  by  excitation 

Figure  4. — Relation  between  vibration  energy  of  gearbox  and 
work  done  by  exciting  force.  (Work  done  is  calculated  from 
the  nonnalized  displacements.) 


Frequency,  Hz 


Figure  5. — Frequency  response  of  gearbox  due  to  unit  exciting 
force.  (Top  plate  is  optimized.) 


operating  at  2700  rpm.  The  initial  vibration  energy  at  this  fre¬ 
quency  is  about  0.976x10'*  J.  The  limits  of  thickness  are  specified 
as  2  mm  ^  t  <  12  mm. 

The  energy  converges  to  a  value  of  about  19  percent  of  the 
initial  value  after  six  iterations  (Fig.  6).  The  frequency  response 
after  optimization  is  shown  by  the  solid  line  in  Fig.  5.  The 
resonance  peak  at  1280  Hz  shifts  to  approximately  1320  Hz.  The 
response  variation  may  look  slight  in  the  figure,  but  the  vibration 
energy  is  reduced.  This  optimum  design  solution  required  about 
875  sec  per  iteration  on  the  Cray  XMP/28. 


Figure  6. — Convergence  of  vibration 
energy.  (Top  plate  is  optimized.) 


Figure  7  illustrates  the  optimum  shape  of  the  top  plate.  The 
thickness  increases  above  bearing  1  and  near  the  edge  that  connects 
to  surfaces  2  and  4.  The  displacement  along  the  Y  axis  at  bearing 
1  is  large  in  the  12th  mode,  1277  Hz  (Fig.  3(b)).  The  increase  in 
thickness  causes  an  increase  of  rotational  rigidity  about  the  X  axis 
to  control  the  displacement.  Figure  8  shows  the  energy  share  of 
each  surface  before  and  after  the  optimization.  The  energy  reduction 
in  the  top  plate  is  comparatively  small.  However,  the  variation  in 
thickness  of  the  top  plate  causes  a  vibration  energy  reduction  in  the 
other  plates. 

Design  of  the  Top  Plate  and  Two  Side  Plates 

In  the  second  example,  two  side  plates  3  and  5  as  well  as  the 
top  plate  are  optimally  designed.  The  side  plates  are  selected 
because  they  do  not  support  the  shafts.  The  conditions  of  excitation 
and  constraint  are  the  same  as  in  the  first  example.  Figure  9  shows 
the  frequency  response.  The  energy  is  considerably  reduced  between 
1180  to  1370  Hz,  and  the  vibration  energy  reduction  is  larger  than  in 
the  first  example.  The  energy  is  reduced  sharply  at  the  first  itera¬ 
tion  and  converges  on  the  value  of  approximately  3  percent  of  the 
initial  value  (Fig.  10). 
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Gearbox  face 

igure  8.— Share  of  vibration  energy.  (Top  plate  is  optimized.) 
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Figure  1 0. — Convergence  of  vibration 
energy.  (Top  plate  and  two  sides 
are  optimized.) 


Frequency,  Hz 


Figure  9.— Frequency  response  of  gearbox  due  to  unit  exciting 
force.  (Top  plate  and  two  sides  are  optimized.) 


The  optimum  shape  of  the  top  plate  and  two  side  plates  are 
illustrated  in  Fig.  II.  The  thickness  distribution  of  the  top  plate  is 
similar  to  the  result  in  the  first  example.  However,  the  weight  of 
the  top  plate  is  reduced  approximately  15  percent  by  the  optimisa¬ 
tion  as  compared  to  the  initial  weight.  The  optimum  shapes  of  the 
two  side  plates  are  similar  to  each  other.  The  thickness  increases  at 
the  bottom  and  sides  of  the  plates.  The  weight  of  plate  S  is  also 
reduced  approximately  7  jjercent  from  the  initial  weight.  However, 
the  weight  of  plate  S  increases  approximately  23  percent  because  of 
the  constraint  condition  of  constant  weight.  Figure  12  shows  the 
energy  share  of  each  surface.  The  energy  is  reduced  in  every  surface, 
and  the  reduction  in  plate  3  is  substantial.  This  may  be  partly  due 
to  the  increase  in  weight. 

Verification  and  Revision  of  the  Proposed  Method 

An  experimental  verification  may  be  required  for  this  kind  of 
design  problem  for  two  main  reas-  is:  to  verify  the  proposed  method 
and  to  estimate  the  error  of  the  result.  Method  verification  can  be 
done  using  computer  simulation  instead  of  experimentation.  After 
obtaining  the  optimum  shape  of  a  beam,  the  authors  showed  that 
the  displacement  amplitude  was  reduced  considerably  following  opti¬ 
misation  by  solving  the  forced  vibration  problem  of  the  beam  (Inoiie, 
Townsend,  and  Coy,  1992).  This  validates  the  proposed  method  it¬ 
self.  Estimating  the  error  of  the  result  depends  on  the  assumptions 
and  simplifications  of  modeling  discussed  in  the  section  on  modeling 
and  eigenvalue  analysis.  A  comparison  with  the  measured  frequency 
is  presented  also  in  this  section. 

A  revision  of  the  method  could  be  completed  later  and  applied 
to  the  design  of  low  vibration  and  lightweight  gearboxes  by  elimina¬ 
ting  the  constraint  of  constant  weight  and  introducing  appropriate 
side  constraints. 
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Figure  1 1 . — Optimum  shape  of  the  top  plate  and  two  side  plates 
of  gearbox. 


Gearbox  face 

Figure  12. — Share  of  vibration  energy.  (Top  plate  and  two  side 
plates  are  optimized.) 


The  method  proposed  to  design  a  low  vibration  gearbox  is  based 
on  a  finite  element  shell  analysis,  modal  analysis,  and  a  structural 
optimization  technique.  The  overall  vibration  energy  is  adopted  as 
the  objective  function  and  minimized  by  varying  the  element  thick¬ 
ness  under  the  constraint  of  constant  weight.  A  computer  program 
was  developed  for  this  purpose. 

The  computer  program  was  applied  to  the  design  problem  of  a 
test  stand  gearbox  used  at  NASA  Lewis.  The  gearbox  was  excited 
by  unit  harmonic  forces  applied  at  the  bearing  positions  in  the 
direction  of  the  line  of  action.  The  work  done  was  estimated  by  the 
product  of  the  bearing  forces  and  the  relative  displacements  obtained 
from  the  eigenvalue  analysis.  The  work  was  approximately  propor¬ 
tional  to  the  total  vibration  energy  of  the  gearbox.  This  suggests 
that  the  gearbox  vibration  magnitude  can  be  estimated  from  the 
relative  work. 

The  top  plate  and  the  two  side  plates  of  the  gearbox  were  opti¬ 
mized.  The  vibration  energy  was  reduced  at  the  excitation  fre¬ 
quency.  The  vibration  reduction  was  obtained  by  variation  of  the 
thickness  distribution  and  the  consequent  shift  of  weight  from  the 
top  plate  to  the  side  plate.  The  contribution  of  each  surface  on  the 
total  vibration  was  determined.  Optimization  of  the  top  plate  only 
had  considerable  effect  in  reducing  the  vibration. 


Figure  13. — Error  of  eigenvalue  analysis  of  simply  supported 
rectangular  plate. 
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APPENDIX  -  ACCURACY  OF  THE  FINITE  ELEMENT 
ANALYSIS 

The  eigenvalue  problem  of  a  simply  supported  rectangular  plate 
is  solved  by  the  developed  program.  The  plate  is  divided  into 
4x4x2,  6x6x2,  and  8x8x2  elements,  respectively.  The  mesh  pat¬ 
terns  are  similar  to  the  pattern  of  the  top  surface  of  the  gearbox 
shown  in  Fig.  2.  Most  of  the  computed  natural  frequencies  are 
smaller  than  the  analytical  values,  and  the  error  depends  on  the 
mode  shape  or  the  mode  order.  The  squared  error  arithmetic  mean 
is  evaluated  and  its  square  root  is  illustrated  in  Fig.  13.  In  the 
lower  30  gearbox  modes,  which  are  used  for  the  modal  analysis,  each 
surface  has  1  and  3  nodal  lines  at  most.  This  mode  corresponds  to 
the  10th  mode  of  the  plate  shown  in  Fig.  13,  suggesting  that  the  FE 
mesh  in  each  surface  of  the  gearbox  may  be  used  for  the  calculation 
of  the  mode  shape  with  an  error  of  approximately  10  percent.  This 
does  not  directly  show  the  error  of  the  gearbox  vibration  because  of 
its  complexity. 
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